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Hard X-ray sources from femtosecond (fs) laser-produced plasmas, including the betatron X-rays from laser 
wakefield- accelerated electrons, have compact sizes, fs pulse duration and fs pump-probe capability, making 
it promising for wide use in material and biological sciences. Currently the main problem with such betatron 
X-ray sources is the limited average flux even with ultra-intense laser pulses. Here, we report ultra-bright 
betatron X-rays can be generated using a clustering gas jet target irradiated with a small size laser, where a 
ten-fold enhancement of the X-ray yield is achieved compared to the results obtained using a gas target. We 
suggest the increased X-ray photon is due to the existence of clusters in the gas, which results in increased 
total electron charge trapped for acceleration and larger wiggling amplitudes during the acceleration. This 
observation opens a route to produce high betatron average flux using small but high repetition rate laser 
facilities for applications. 

Synchrotron light sources have proven their usefulness for users especially in the biological and condensed 
matter applications 1 . Moreover, some X-ray Free Electron Laser (XFEL) facilities are recently constructed 
and operational 2 , which are expected to have revolutionary impacts on science and technology. However, 
these machines are huge in size (km-size) and accessible to limited users. Recently with the development of 
ultrashort high power lasers, the laser-plasma X-rays are attracting increasing interest as compact and affordable 
sources with some unique features for many applications. Hard X-ray emission from femtosecond laser-produced 
plasmas has been extensively studied in the past years 311 . The source energy conversion efficiency and temporal 
duration have been greatly improved 6 ' 7 . However, most laser-driven hard X-ray sources are spatially symmetric, 
except for high harmonic generation 8 , Thomson scattering 9 and betatron radiation 3 . High harmonics and 
Thomson scattering usually suffer from their low yield. 

A well-collimated X-ray beam with fs duration is produced when a relativistic laser pulse propagates in an 
under dense plasma: this occurs as some electrons are injected into the laser wakefield in plasma and experience 
transverse betatron oscillations while being accelerated forward 3 . Intense betatron radiation, with photon energy 
from keV X-rays 10 to gamma-rays 11 , has been generated via the accelerated electrons undulating in the bubble 
structure of the wakefield. These oscillations occur at the betatron frequency cop = w p /(2y) 1/2 , where w p is the 
plasma frequency and y is the Lorentz factor of the electron beam 10 . The average photon number with mean 
energy fjco c emitted by an electron is given by Nx = 5.6 X 10~ 3 NpK, where Np is the number of oscillation 
periods and K is the strength parameter of the plasma wiggler given by K = 2n(yr 0 )/X b = 1.33 X 
10" 10 y 1/2 n e 1/2 [cm~ 3 ]r 0 [inm], in which r 0 is the electron wiggling amplitude. Hence the characteristics of betatron 
X-rays are mainly determined by the accelerating electron properties, including the electron energy, beam charge, 
as well as the electron wiggling amplitude inside the wakefield structure. However, when using a gas target, it is 
usually mutually exclusive to improve both the accelerated electron energy and charge simultaneously. For 
example, a low plasma density is necessary for longer pump -depletion length and higher electron energy, but 
it also results in lower total electron charge being injected. This is the reason why the reported electron charge in 
most of experiments using a medium-size laser facility is limited to 100 pC. To enhance the betatron X-rays, an 
experiment attempted to increase the wiggling amplitude of injected electrons using an asymmetric tilted pulse 
front 12 , but it usually effects the formation of the laser wakefield and reduces the final electron energy. Therefore, 
in order to enhance the betatron radiation brightness, one needs to find an effective way to increase both the 
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Figure 1 | Schematic of the experimental setup. 



electron beam charge and the electron wiggling amplitude in laser 
acceleration, which has been the subject of many current investi- 
gations. 

In this paper, we present the generation of large-charge electron 
beams as well as strong betatron X-rays using an Ar clustering gas jet 
target irradiated with a 3 TW femtosecond laser. The measured 
emission yield with the photon energy > 2.4 keV reaches 2 X 10 8 
photons/shot, which is ten-fold enhancement comparing to the 
emission flux produced by using gas target with the same laser para- 
meters. Experiment and simulations point to the existence of cluster- 
ing in the background gas, which leads to the increased number of 
trapped electrons as well as much larger wiggling amplitudes during 
acceleration, beneficial to betatron X-ray photon emission. 

Results 

The experimental set-up is shown in Fig. 1 . Laser pulses from a 3 TW 
laser system at 80 fs are focused to the target center either with Ar 
clustering gas target or He gas target with the same plasma density 
(see methods). By using Ar clustering gas as an acceleration media, 
the X-ray photon yield of laser-driven betatron radiation is greatly 
enhanced in our experiment. Fig. 2(a) shows the Ar betatron X-ray 
profile obtained in a single laser shot with the photon energy 
> 2.4 keV. Compared to the X-ray profile from a gas He target with 
the same laser parameters, the Ar betatron emission is found to be 
much more intense. The emission yield reaches 2 X 10 8 with a 
confidence interval of 1.6 ~ 2.4 X 10 8 for a single laser shot with a 
divergence angle —10 mrad, whereas it is < 1 X 10 7 in the case of the 
He gas target, as shown in Fig. 2(b). This is an achievement of intense 
keV betatron X-rays using laser system at a few TW only. Fig. 2(c) 
shows the measured X-ray signal for different cut-off filters in experi- 
ment. Most of betatron X-ray photons are concentrated in the energy 
higher than 2.4 keV. Compared to a previous observation 13 with a 
gas target and a 10 TW laser facility, the betatron X-ray flux at the 
similar X-ray photon energy (~ 3 keV) we obtained here is 400 times 
larger even though our laser power is 3 times lower, i.e. ~ 1000 times 
higher conversion efficiency is achieved compared to the previous 
reported results. Using a synchrotron asymptotic limit distribution 
model as defined in Ref. [ 14] , the critical energy we obtained is E crit = 
hco c = 0.6 keV after minimizing the least squares fit for data obtained 
from experiment and the synchrotron model. It should be empha- 
sized that the betatron X-ray yield depends critically on the laser 
contrast. The yield is reduced by two orders of magnitude when 
the laser pulse contrast decreases from 10 9 to 10 7 while fixing the 
laser pulse energy. Pre-expansion of a solid-density cluster by the 
laser pre-pulse leads to ineffective electron heating and acceleration, 
as described below, and results in the decrease of X-ray yield 
produced. 

The accelerated electron spectrum was measured simultaneously 
because it reflects the behavior of the electron wiggling that leads to 
betatron X-ray generation. As shown in Figs. 3(a) and 3(b), the 
electron beam profiles are measured using an Al foil -wrapped 
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Figure 2 | The measured betatron emission characteristics: The emission 
beam profile obtained using 36 um Al foil wrapped IP for Ar target (a) and 
He target (b), respectively. The hallow cloud signal in (a) (left side) was 
caused by low energy electrons bypassing the magnet. The measured 
(square) and modeled (triangle) X-ray signal for different filters are shown 
in (c). The inset shows the beam profile recorded after the cut-off filters 
18 umAl(I), 18 umAl+ 10 umCu(II), 18 umAl + 20 um Cu (III), and 
18 um Al + 30 um Cu (IV), with different cutoff energy 2.4 keV, 3.4 keV, 
4.1 keV and 6.2 keV, respectively. Modeling is based on a best fit of 
synchrotron spectrum with E crit = 0.6 keV. 
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Figure 3 | Electron beam characteristics: The beam profile obtained using 
18 Jim Al foil wrapped IP for He target (a) and Ar target (b), respectively, 
with the same electron density (2 X 10 19 cm" 3 ) and the same laser 
parameters. The electron spectra are shown in (c) and (d) obtained after 
the dispersive magnet. The color scale is the same in (c) and (d). 
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imaging plate directly after the nozzle (without magnet) in the cases 
of an Ar clustering gas target and a pure He gas target. The beam 
charge reaches over 200 pC with electron energy higher than 1 MeV. 
It is over 50 times higher than the case of He gas with the same 
electron density used (2 X 10 19 cm" 3 ). We note that the divergence 
angle of the electron beam using Ar cluster target is — 17 mrad, 
much larger than ~ 1 mrad for the case of He gas target. It indicates 
that the electron is not well collimated after acceleration in the clus- 
tering case, which suggests that large betatron oscillations of electron 
may occur during the electron acceleration. In addition, because of 
the large charge observed with the Ar cluster target, the space charge 
effect may play some role in the early stage of electron acceleration. 
Figs. 3(c) and 3(d) show the spectra obtained with the same inter- 
action conditions. A quasi- monoenergetic electron beam at the cent- 
ral energy about 21 MeV is obtained for the He-gas case. But in 
contrast, for the case of Ar cluster target, we observe a continuous 
spectrum with the maximum energy over the 20 MeV level. 
Furthermore, both the beam charge and beam divergence angle in 
the case of Ar cluster are much larger than the case of He gas target. 
This measurement is consistent with the beam profile results. It 
shows the clustering gas target is more suitable for obtaining high- 
charge electron beam acceleration, than the gas target, with larger 
divergence angle. We use these electron parameters to estimate the 
betatron X-ray radiation. The critical energy E crit [keV] = 5.3 X 
10" 24 y 2 n e [cm" 3 ]r 0 [|im] = 0.53 keV, which is similar to our experi- 
mental fitting above. The plasma wiggler strength K = yk b r 0 =14.1 
for photons with energy > 2.4 keV. For a betatron wavelength ^ b = 
(2y) 1/2 ^ p = 44 urn, the number of betatron oscillations executed by 
an electron for this laser-produced ion channel N 0 = cT int /X h = 22, 
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Figure 4 | Simulation results. 2D snapshots of the electron density 
distribution of gas (a) and clusters (c) at time of t = 140T 0 (T 0 is the laser 
optical cycle), (e) shows the electrical field in the interaction regime at the 
same time [corresponding to (c)] for the cluster target. The spatial 
distributions of the accelerated electrons at the same time are shown for the 
case with the gas target with energy > 10 MeV (b) and clusters with energy 
> 10 MeV (d) as well as > 20 MeV (f). (g) shows the electron energy 
spectra obtained with gas (blue line) and clusters (red line). 



where cT int is the interaction length which should be about 1 mm in 
the experiment. So, the total betatron yield with energy of 2.5 keV is: 
N x « 4.4 X 10- 4 N e N 0 K = 1.7 X 10 8 . This value is in reasonable 
agreement with experimental detection. 

Discussions 

To further clarify the involved physics, two-dimensional (2D) fully 
electromagnetic particle-in-cell (PIC) simulations have been per- 
formed. Figures 4(a) and 4(b) show snapshots of the density distri- 
bution of total electrons and energetic electrons (> 10 MeV), 
respectively, at t = 0.43 ps for the case of a pure gas target. Small 
amounts of electrons are injected into the second and following 
wakefield wave buckets and the tightly transverse collimation of 
the electron beam indicates that the trapped electrons experience 
small oscillation amplitudes during acceleration, which is regularly 
observed in self-modulated laser wakefield acceleration (SM- 
LWFA) 15 . In this case, several quasi-monoenergetic electron bunches 
with different energies can be achieved, also shown in Fig. 4(g) 
(simulation) and Fig. 3(c) (experiment) of the electron spectra. 
However, when using a clustering gas target, the electron density 
distribution shows a large bow-like structure, as shown in Fig. 4(c) 
for the snapshot at t = 0.43 ps, and the electron density exhibits 
serious modulation in a broad region directly following the laser 
pulse. A big size plasma cavity 16 (~ 10 X 0 ) appears in the center area, 
which heats the most energetic electrons, as seen in Fig. 4(e). A large 
quantity of electrons with the charge of 300 pC is accelerated and the 
electron beam transverse profile becomes very broad, as shown in 
Figs. 4(d) and 4(f). The total accelerated electron charge in the case of 
clustering gas target is seen to reach ~ 170 pc for electrons with 
energies > 20 MeV. However, in case of a pure gas, electron charge 
injected into wakefields and accelerated to the same energy is much 
lower. Note that during the propagation of the laser pulse in the 
clustering gas, the laser pulse develops a wing structure, which also 
leads to a wing structure in the accelerated electron beam. 

Our simulation results show that the direct laser acceleration 
(DLA) 17,18 has important contribution to drive electrons with large 
oscillations in the case of a clustering gas target. Figure 5(a) plots 
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Figure 5 | (a) Traced electron trajectories from simulation with a cluster 
target (solid line) showing transverse deflection y. Two dotted lines for the 
case of pure gas target are also shown for comparison, (b) A traced electron 
average energy versus time, where the solid red line shows the energy 
gained from DLA and the dashed blue line shows that from the ES field, (c) 
Energy gain from DLA versus that from ES field at 200 T 0 , where each point 
denotes a traced electron; here black dots are found with the Ar clustering 
target and red dots with the He gas target. 



SCIENTIFIC REPORTS | 3 : 1912 | DOI: 1 0.1 038/srep01 91 2 



3 



some trajectories of selected electrons which are accelerated to high 
energy. It shows that the oscillation amplitudes increase from a cen- 
ter cavity quickly to 10 um and electrons undergo wiggling at the 
betatron frequency (Op, indicating a strong resonant coupling with 
the laser field. As a comparison, the electron transverse oscillation is 
very small in the case of He. In the meantime, the transverse 
momentum will transfer to longitudinal momentum via the vx B 
drift by a self-generated azimuthal magnetic field, where electrons 
gain energy directly from the laser. For example, in Fig. 5(b), the 
energy gain of a sample electron reaches — 30 MeV over several 
betatron oscillations. It also indicates that the energy obtained from 
transverse electrical field is more important than that achieved from 
the longitudinal electrical field, also see in Fig. 5(c). We count the 
total energy gain for all electrons with energy over 10 MeV. It turns 
out that the energy gain from DLA counts for 76% of the total energy 
gain. In contrast, in the case with the He pure gas target, the energy 
gain is mostly due to the laser wakefield acceleration (only about 
1.35% of the total energy gain comes from DLA). The energy spec- 
trum is in thermal distribution, as Fig. 4(g) shown, because of the 
random original acceleration phase of accelerated electrons. On 
the other hand, the electron oscillation amplitude as simulated 
(—10 |im) is in good agreement with that deduced from the critical 
energy (0.6 keV) fit to the experimental betatron X-ray spectrum, — 

9 urn. All of the above simulation and experimental results strongly 
suggest the presence of DLA mechanism in the case of clustering gas 
target. Usually when using pure gas target, DLA will dominant over 
LWFA at high plasma density and very high laser intensity, i.e: when 
a 0 = 0.89 X I 1/2 [10 18 W cm- 2 ]^ 0 [|im] > 1, which always take place 
by using petawatt laser pulses 14 . However, as described in Ref. [19], 
when using a clustering gas target even in modest laser intensity, the 
strong resonant heating of cluster electrons 7 in the early stage of 
interaction results in higher local electron density and much larger 
pre-acceleration phase in the transverse direction than the case of a 
non- clustering gas, which makes it easy to shorten the time period 
before strong vxB starts. As shown in Fig. 5(a), electrons start to be 
heated from the very beginning of the interaction and quickly experi- 
ence betatron resonant oscillations. Therefore, the clustering gas is 
an ideal medium to stimulate the DLA mechanism under modest 
laser intensity. 

In summary, we have demonstrated a new method for generat- 
ing intense betatron X-rays using a clustering gas target irradiated 
with an ultra-high contrast laser of 3 TW only. The yield of the 
well-collimated Ar X-ray betatron emission has been measured to 
be 2 X 10 8 photons/pulse, which is a ten-fold enhancement com- 
pared to the emission flux produced by using a normal gas target 
under the same laser parameters. Simulations point to the exist- 
ence of clustering as a contributor to the DLA mechanism, lead- 
ing to higher accelerated electron charge (x50) and much larger 
electron wiggling (—10 urn) amplitudes in the plasma channel, 
thereby finally enhancing the betatron X-ray photons. Here in 
this DLA regime, electrons are efficiently driven on the fs time 
scale, producing enhanced X-ray emission with durations about 
three orders of magnitude shorter than that of typical pulses 
produced by synchrotron sources. Together with the inferred 
source size of the betatron oscillation amplitude —10 um and 

10 fs duration deduced from electron bunch length in Fig. 4(f), 
the peak brightness of the radiation in the energy range of about 
3 keV is estimated to be — 5 X 10 21 photons/s/mm 2 /mrad 2 / 
0.1 %BW, which is comparable to the peak brightness of the third 
generation of synchrotron radiation sources. Being perfectly syn- 
chronized with the driver laser pulse, this kind of table-top, ultra- 
intense and spatially coherent X-ray sources with fs duration 
opens prospect for time-resolved "single-shot" ultrafast pump- 
probe applications employing phase contrast imaging, X-ray 
absorption fine structure and X-ray diffraction without heat- 
induced blurring. Moreover, a variety of other applications 



including imaging with high spatial resolutions may be realized 
with such X-ray sources at much higher photon flux produced 
with 0.1-1 kHz lasers. 

Methods 

Laser system. The experiments are carried out by using the XL- II laser facility in the 
Institute of Physics -CAS, which includes a 10 Hz, 300 mj Ti:Sapphire laser working at 
the centre wavelength X = 800 nm. The pulse with duration x 0 = 80 fs is focused with 
an f/6 off-axis parabola (OAP) onto a focal spot of size w 0 = 4 ± 0.1 um. The focal 
spot size was measured by a removable diagnostic system consists of a microlens and a 
visible CCD. In the focal region the laser average intensity is I = 3.0 X 10 18 W/cm 2 . 
With the help of optical parametric chirped pulse amplification, the laser pulse 
contrast, compared to its ns prepulse, has been improved to 10 9 . 

Production of atomic clusters. A supersonic pulsed gas (Ar) jet is used, which is a 
1.2 mm X 10 mm rectangular nozzle at the exit. With sufficiently high backing 
pressure, clusters are formed in the gas jet flow due to the adiabatic cooling of a gas 
expanding into vacuum. A probe beam is used for detecting the parameters of the 
clusters. According to the Hagena scaling law 20 , the cluster size is estimated. An 
average size of ~ 8 nm in diameter is estimated at a stagnation pressure of 4 Mpa, 
with an average plasma density of 2 X 10 19 cm -3 . 

Diagnostics of electron and radiation. Electrons are dispersed by a strong magnet 
placed 20 cm after the nozzle, and an imaging plate follows the magnet to record the 
electron beam. A filtered imaging plate is also located in the laser propagation 
direction to detect the X-ray flux and spatial profile with photon energy > 2.4 keV. 
Different filters will be placed in front of the IP (as described in Fig. 2) to selectively 
attenuate the X-ray flux and estimate the spectral content as cut-off filters 14 . 

Particle-in-cell simulation. Numerical simulations were carried out with our PIC 
code KLAP 2D, where a p-polarized laser pulse is focused into a pure gas target or a 
cluster-gas target with a focal spot size w 0 = 8 um. The laser pulse has a transverse 
Gaussian shaped envelope with a peak strength a 0 = 2.5, corresponding to an average 
intensity I = 3 X 10 18 W/cm 2 in the focal region. The plasma density used in 
simulation for He is 2 X 10 19 cm -3 and for Ar surrounding gas is 1.7 X 10 19 cm -3 . The 
clusters, with a fixed density n = 1.7 X 10 22 cm -3 and diameter d = 8 nm, are 
distributed uniformly in the background gas with a fixed gap r c = 120 nm between 
one another. Simulations were performed in a boosted frame moving along with the 
laser pulse. The size of the simulation box is 36 X 80^ 0 2 , with grid size of 0.0 1^ 0 
(9 particles per cell). To verify that electrons gain energy mainly from the laser fields 
rather than from the longitudinal fields, one chose some particles with final energy 
E > 30 MeV to trace their energy gains. According to the equation of motion for 
electrons, for each particle the energy gain from electrostatic field and laser field 
(DLA) can be obtained from E ES = Jo'ExVxdt' and E DLA = Jo'ExV^dt' , respectively, 
where E x and E ± are the normalized longitudinal and transverse electric fields, 
respectively 21 . 
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